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Predicting the magnetic vectors within coronal mass ejections 
arriving at Earth: 1. Initial Architecture 
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Thompson,^ A. Pulkkinen,^ R. Evans,^ T. Nieves-Chinchilla,^’^ 

Abstract. The process by which the Sun affects the terrestrial environment on short 
timescales is predominately driven by the amount of magnetic reconnection between the 
solar wind and Earth’s magnetosphere. Reconnection occurs most efficiently when the 
solar wind magnetic field has a southward component. The most severe impacts are dur¬ 
ing the arrival of a coronal mass ejection (CME) when the magnetosphere is both com¬ 
pressed and magnetically connected to the heliospheric environment. Unfortunately, fore¬ 
casting magnetic vectors within coronal mass ejections remains elusive. Here we report 
how, by combining a statistically robust helicity rule for a CME’s solar origin with a sim¬ 
plified fiux rope topology the magnetic vectors within the Earth-directed segment of a 
CME can be predicted. In order to test the validity of this proof-of-concept architecture 
for estimating the magnetic vectors within CMEs, a total of eight CME events (between 
2010 and 2014) have been investigated. With a focus on the large false alarm of Jan¬ 
uary 2014, this work highlights the importance of including the early evolutionary ef¬ 
fects of a CME for forecasting purposes. The angular rotation in the predicted magnetic 
field closely follows the broad rotational structure seen within the in situ data. This time- 
varying field estimate is implemented into a process to quantitatively predict a time-varying 
Kp index that is described in detail in paper II. Euture statistical work, quantifying the 
uncertainties in this process, may improve the more heuristic approach used by early fore¬ 
casting systems. 


1. Introduction 

CMEs are often observed to have twisted “flux rope” mag¬ 
netic field structures [Liu et al.^ 2008; Vourlidas^ 2014]. If 
favorably oriented, these can lead to extended southward 
excursions of the interplanetary magnetic field (IMF) as 
the CME passes by, resulting in periods of enhanced re¬ 
connection on Earth’s dayside and energy input into the 
magnetosphere. Draping of the IMF around the CME as it 
moves through the solar wind may also give rise to south¬ 
ward fields. In contrast, northward-directed fields inhibit 
reconnection, resulting in a weaker magnetospheric response 
[Dungey^ 1961]. While prolonged southward fields are often 
observed without the presence of a clear structured tran¬ 
sient, the additional plasma parameters often associated 
with CMEs usually make them the most geo-effective events 
[e.g., Tsurutani and Gonzalez, 1997; Zhang and Moldwin, 
2014]. Thus, inferring the direction of the flux rope fields 
inside a CME before it arrives at Earth would be a major 
advance in geomagnetic activity prediction. 

In addition, the CME’s initial configuration and its in¬ 
teraction with the inhomogeneous ambient solar wind can 
lead to deformations, rotations, and deflections of the mag¬ 
netic field, which are difficult to quantify [e.g., Odstrcil and 
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Pizzo, 1999; Savani et ah, 2010; Nieves-Chinehilla et ah, 
2012]. Distortions of CMEs have previously been observed 
by coronagraphs. However, their influence on the magnetic 
structure is difficult to estimate because the magnetically- 
dominated regions of CMEs appear as dark cavities within 
images, such as those seen by the STEREO spacecraft 
[Howard and DeForest, 2012]. Therefore, a common ap¬ 
proach to predicting magnetic vectors within a CME prop¬ 
agating towards Earth is to use solar observations as inputs 
into 3D computational simulations. Unfortunately, obtain¬ 
ing realistic magnetic field directions at Earth from such 
calculations is scientifically challenging and computationally 
intensive [Manehester et ah, 2014]. 

Thus, models used for routine CME forecasts by various 
space weather services do not include magnetic structures 
within the simulated CMEs [e.g., Zheng et ah, 2013; Shiota 
et ah, 2014]. For example, ENLIL models the propagation 
of CMEs from ~ 20 solar radii (Rs) to beyond Earth at 215 
Rs and includes the background solar wind magnetic field. 
However, the CME is simplified to a high pressure plasma 
pulse with a size and propagation direction estimated from 
solar imagery [Zheng et ah, 2013]. CME arrival-time predic¬ 
tions from these models provide lead times of ~ 2 — 3 days, 
and their accuracy has been well investigated [Taktakishvili 
et ah, 2010; Vrsnak et ah, 2014; Colaninno et ah, 2013]. In 
contrast, the important magnetic vector information is only 
revealed when in situ measurements are made by spacecraft 
upstream of Earth at the first Lagrangian position (LI) ~ 1 
hour prior to the CME arriving at Earth, thereby severely 
limiting the lead time available for reliable, magnetic field- 
based, storm warnings. 

Difficulties in observationally determining the magnetic 
profile of a CME arriving at Earth from only solar imagery 
predominately lie with several complex stages that change 
the initial solar configuration to the final topological struc¬ 
ture at Earth. We suggest that for forecasting purposes, 
statistically significant predictions can be made by simplify¬ 
ing the complex behavior to a core set of parameters. In this 
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paper, we highlight three key components of a proof of con¬ 
cept developed to improve the prediction of a storm’s sever¬ 
ity: 1. the use of the hemispheric helicity rule to provide a 
robust initial magnetic configuration at the Sun; 2. Define 
a ‘volume of influence’ of the CME, within the heliosphere, 
for which the Earth’s trajectory can be estimated; and 3. 
incorporating magnetic vectors from a simplified magnetic 
flux rope model to create a time-series upstream of Earth. 

Erom the analysis of eight Earth-directed CMEs between 
2010 and 2014, we conclude that the incorporation of mag¬ 
netic field vectors in this way can significantly improve ge¬ 
omagnetic forecasts by providing a time-varying magnetic 
profile of the CME. The time varying magnetic profile is 
then incorporated with an experimental technique to create 
a time varying Kp index forecast that replicates the forecast 
deliverables by NOAA. Eurther details on the geomagnetic 
indices and their uncertainties that are borne from the vec¬ 
tor estimates are described by Savani et al. (2015, herein 
referred to as Paper II). 

2. Event selection 

This article discusses the proof-of-concept architecture for 
estimating the magnetic vectors with the aid of a case study 
CME event that was released from the Sun on January 7th 
2014 (see Eigure 1). The recent release of this event allows 
for comparisons between the results described here and the 
current processes employed by real time space weather fore¬ 
casters and their estimated geomagnetic indices (described 
further in Paper II). 

A total of eight Earth-directed CME events between 2010 
and 2014 were selected with three driving criteria: 1. Un¬ 
ambiguously define the solar source of the overlying field 
arcade from a single or double active region and possibly 
with an eruptive flare (see section 3 for more details); 2. A 
clear leading edge structure from multiple remote observa¬ 
tions to unambiguously define the size and orientation of the 
CME morphology (see section 4 for more details); and 3. a 
significant measurable effect by geomagnetic indices. 

The eight events described in this paper were chosen from 
a CME list compiled by Colaninno et al. [2013] and Pat- 
sourakos, S.(personal list), with details of each event dis¬ 
played in Table 1. Eurther Earth-directed CME lists with 
more generic requirements have also been published [e.g., 
Richardson and Cane, 2010; Mostl et al, 2014]. 

The hemispheric solar source region of the CME is iden¬ 
tified from solar observations. Eigure 1 displays a 17lA 
image from the AIA instrument onboard the SDO space¬ 
craft [Lemen et al, 2012] taken at 20.14 UT on 7th Jan¬ 
uary, 2014. This event has an inconclusive Earth-arrival 
time and in situ profile, and has been chosen to highlight 
the complexity in forecasting processes. The uncertainty 
from this event stems from the predicted arrival time being 
approximately 24 hours earlier than when on-duty forecast¬ 
ers labeled the actual arrival from real-time LI in situ data. 
Throughout this period, the solar wind plasma parameters 
displayed significantly lower velocities than were expected 
as well as missing a strong and distinctive magnetic field 
rotation of an obstacle. 

3. Solar Initiation 

The helicity and initial orientation of the magnetic 
flux rope structure within a CME are inferred from the 
“Bothmer-Schwenn” scheme. This relates the flux rope 
properties to sunspots, the solar cycle, and whether the 
CME originates on the northern or southern solar hemi¬ 
sphere [Bothmer and Schwenn, 1998]. The reliability of 
this solar hemispheric rule remains controversial. It was 


only in late-2013 when the probability of a CME’s topology 
conforming to the hemispheric rule was re-confirmed to be 
>80% [Wang, 2013; Hale, 1925]. Thus, the initial helicity 
and field structure of CMEs can be inferred from this scheme 
with a reliability that is likely to be ~80%. 

Ordinarily, a CME is linked to a single active region where 
the standard Bothmer-Schwenn scheme should be applied. 
However in cases such as this January 2014 event, the mag¬ 
netic loop structure before eruption has a leading negative 
polarity spanning over two active regions. Thus, a South- 
West-North, “SWN”, flux rope field direction from southern 
hemisphere of solar cycle 23 under the Bothmer-Schwenn 
scheme is appropriate. This implies the CME has a right 
handed chirality. 

Harra et al. [2007] highlighted the complexity of estimat¬ 
ing the orientation of an interplanetary CME from simple 
solar observations. The work displayed that two CMEs re¬ 
leased in November 2004 from a similar source location had 
drastically different final topologies. However, this can be 
reconciled with the Bothmer-Schwenn scheme if the differ¬ 
ent polarity of the active region’s leading edge is taken into 
account. 

In this article, we consider six simpler CMEs released 
from a single active region and examine whether it is possible 
to generate more reliable predictions of the field structure at 
1 AU. We also investigate two more complicated cases where 
connected active regions are involved (September 27th, 2012 
and January 7th 2014). 

4. Remotely Sensed Evolution 

Since deflections, rotations and other interactions may oc¬ 
cur during CME propagation to Earth, the initial Bothmer- 
Schwenn configuration is adjusted using coronagraphic data 
from the SOHO and STEREO missions [Brueckner et al, 
1995; Howard et al, 2008]. The final tilt and source region 
of the magnetic flux rope, after which radial propagation 
is assumed [Nieves-Chinchilla et al, 2013], is estimated us¬ 
ing the graduated cylindrical shell (CCS) model [Thernisien 
et al, 2009], when the CME reached ~ 15 Rs. Eigure 2 dis¬ 
plays images from the COR2 instruments onboard the two 
STEREO spacecraft (A and B) and the LASCO instrument 
onboard SOHO that are used to triangulate the CME struc¬ 
ture. Where three well-separated observations exist, the 
CCS model provides relatively well-constrained estimates of 
the orientation and size of the CME without any ambiguity 
[Liu et al, 2010a; Rodriguez et al, 2011]. The CCS model 
may still be implemented without multi-point observations 
in the same way as various other cone-structure methods 
can be implemented. However in such cases, the possible 
degeneracy in the observational morphology limits all meth¬ 
ods and thus highlights the difficulties in performing reliable 
forecasts. 

The outputs from the CCS model along with estimates 
of the average CME size [Yashiro et al, 2004] are used to 
create a “volume of influence” defined as the volume the 
CME is expected to traverse as it propagates through the 
heliosphere. The shaded region in Eigure 1 displays the pro¬ 
jection of this “volume of influence” onto the Sun, suggesting 
that the Earth grazed the northern edge of this case study 
event. The projected area is calculated from the ‘shadow’ of 
the CME that is assumed to be cylindrical in shape with cir¬ 
cular cross-section. Two parameters (flux rope axis length 
and flux rope width) are required to estimate the projected 
area. The axis length, shown as a dashed curve on Eigure 1, 
is estimated from the half angle, auaw^ given by half the an¬ 
gular width of the CME in a direction parallel to the CCS 
model axis. The projected width of the CME transverse 
to this axis is assumed equal to the average CME width, 
as found from statistical studies [Nieves-Chinchilla et al, 
2013]. 
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Any uncertainty in the inferred CME orientation is likely 
to have only a minimal effect on the predicted magnetic field 
vectors since it is likely to be eclipsed by the larger uncer¬ 
tainties arising from estimating the magnetic field strength 
and the assumption of a symmetric cylindrical flux rope, as 
explained below. Eurther testing of these assumptions are 
addressed in paper II and are tested relative to predicted 
estimates of Kp. The coronagraphic images show how the 
coronal magnetic loops seen in SDO have been deflected to 
the south west (Eigure 1). Here, we use coronagraphic im¬ 
agery to estimate the final CME radial trajectory but future 
work could attempt to increase prediction lead times by, for 
example, incorporating CME deflections by coronal holes 
[Cremades and Bothmer, 2004; Mdkeld et a/., 2013]. 

The shortest (perpendicular) distance between the 
Earth’s projected location and the flux rope axis is indi¬ 
cated on Eigure 1 with a blue curve. Normalizing this to 
the total perpendicular distance to the flux rope outer edge 
(flux rope radius) gives a quantity that is correlated with 
the impact parameter (YO) which is a key parameter for 
in situ flux rope modelers. The theoretical model of the 
impact parameter used in this study is displayed in Eigure 
3; the Earth’s projected arc distance is displayed in solar 
radii. This theoretical function is justified by using a simple 
linear correlation between the Earth’s distance and YO for 
the inner core region surrounding the flux rope axis (inner 
highlighted area in Eigure 1). The outer area is correlated 
with a trigonometric sine function and is designed to physi¬ 
cally represent the distortions to the idealized flux rope that 
occur during propagation as well as possible draping of the 
surrounding solar wind magnetic field outside the actual flux 
rope structure. These distortions to the flux rope are some¬ 
times termed ‘pancaking’ [Riley and Crooker, 2004; Savani 
et a/., 2011a] with recent studies suggesting the inner core 
of a CME is likely to maintain a quasi-cylindrical structure 
while the outer structure may become severely deformed by 
the ambient medium [Demoulin and Dasso, 2009; Savani 
et al, 2013a]. 

5. In Situ Flux Rope 

To generate an estimate time-series of the magnetic vec¬ 
tor direction passing over a fixed point such as LI, we must 
employ a methodology to create a 1-D (spacecraft) trajec¬ 
tory through a theoretical structure, and to define the start 
time of the object at this fixed point. 

5.1. Time of Arrival 

To improve the time of arrival prediction of a CME is 
beyond the scope of this work, and several advances on this 
topic have been performed. Currently there are several pro¬ 
cedures to calculate the speeds of remotely-observed CMEs, 
quantify their deceleration, and forecast their speeds up¬ 
stream of Earth at LI [see further literature within, e.g., 
Owens and Cargill^ 2004; Colaninno et aL, 2013; Tueker- 
Hood et aL, 2015]. We choose to assume a simple average of 
the measured CME speed close to the Sun as determined by 
the NOAA Space Weather Prediction Center (SWPC) and 
the predicted speed at Earth. In the case of the 7th January, 
2014 event, this gives a CME speed of 1300 km/s and a pre¬ 
dicted arrival time of 8th January, 21.45 UT. By combining 
this information with the flux rope model described below, 
a time-series of magnetic vectors is created. 

In order to compare the accuracy of the modeled mag¬ 
netic vector time-series with data and test the technique in 
the ‘research domain’, we manually adjust the arrival time 
of the model fit to the best guess estimate within the LI 
data. The field rotations between the model estimate and 
data were then manually inspected. 

However, for a readily implement able process for estimat¬ 
ing the magnetic vectors in advance, different forecasting 


systems can simply employ their best estimate of the arrival 
time. 

5.2. Flux Rope Model 

The configuration of the magnetic flux rope is calculated 
by assuming a constant-alpha force-free (CAEE) flux rope 
[Burlaga, 1988; Zurbuehen and Riehardson, 2006, and refer¬ 
ences therein] and a cylindrical geometry locally around the 
Earth’s predicted trajectory through the CME. Previously, 
triangulation of the CME direction from remote sensing have 
provided adequate information as to the expected structure 
arriving at LI [Liu et a/., 2010b]. However a Grad-Shafranov 
reconstruction technique used by Liu et al [2010b] would not 
be appropriate in creating a model to estimate the struc¬ 
ture. Euture work may consider implementing a more com¬ 
plex model that better represents the distortions occurring 
to a CME at LI [e.g., Maruhashi, 1997; Hidalgo et al, 2002; 
Owens et al, 2006] 

The magnetic vectors generated along the Earth’s tra¬ 
jectory from a CAEE flux rope model is created from the 
MHD momentum equation under magnetostatic equilib¬ 
rium; which can be reduced to j = aB. A solution to this 
equation can be used to generate a cylindrical magnetic flux 
rope with circular cross section, with the components of the 
magnetic field vector expressed by Bessel functions, and a 
commonly set to 2.41 [Savani et al, 2013b]. Euture work 
can consider reducing a as a simple solution to potential 
flux erosion occurring to the CME during propagation [Ruf- 
fenaeh et al, 2012]. 

The projected axis onto a 2D plane of the CME is pro¬ 
vided by a single angle orientation (0) estimated from CCS 
model. However, the component of the flux rope axis par¬ 
allel to the radial direction is estimated theoretically, by 
measuring the shortest distance between the Earth trajec¬ 
tory through the CME away from the CME nose, Ln. In 
practice, this was performed by measuring half the length 
of the flux rope axis (Rax) and the length between the flux 
rope axis center and the position where the Earth perpen¬ 
dicular position (Eigure 1, blue curve) meets the flux rope 
axis (De); thereby defining, Ln = DeIRux- Ln = 0 repre¬ 
sents the case where the CME nose is propagating directly 
towards Earth, and there is no radial contribution to the 
flux rope axis. However when Ln = 1, the flux rope axis is 
entirely radial in direction, as might be the case when the 
Earth’s trajectory is along a CME leg. Eigure 4 displays how 
the radial contribution to the axis vector is estimated from 
an angular value (A) that varies between 0° (CME nose) 
and 90° (CME leg) in a scheme similar to that expressed by 
Janvier et al [2013]. Both 0 and A are used to create a 3D 
flux rope axis direction. 

The magnitude of the magnetic field along the central flux 
rope axis is assumed in this case study to be 18.0 nT. This 
is calculated by assuming the maximum estimated magnetic 
field strength within the plasma pile-up region simulated by 
the WSA-ENLIL+Cone model (10.3 nT) corresponds to the 
magnetic field strength at closest approach within the flux 
rope structure. The impact parameter obtained using the 
‘volume of influence’ (set at 0.91 for the January 2014 event) 
is then used to estimate the maximum field strength along 
the central flux rope axis. In effect, this technique estimates 
the I B I of a CME from a correlation of the inner he¬ 
liospheric CME velocity and a simulated background solar 
wind field strength driven by magnetograms. 

The flux rope axis direction, chirality, magnetic field mag¬ 
nitude and impact parameter provide a complete set of pa¬ 
rameters to generate a time series of magnetic vectors along 
a theoretical Earth trajectory (Eigure 5 red curves). 
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5.3. Magnetic Field Strength 

The field strength is inferred from a model currently used 
for forecasting by CCMC, so this method could be imple¬ 
mented using existing forecasting capabilities. In the future, 
other methods whose uncertainties have not yet been statis¬ 
tically quantified might be used, for example estimating the 
poloidal and total flux content of a CME from flare rib¬ 
bon brightening [Longcope et al, 2007]; flux rope speed and 
poloidal flux injection to estimate field strength [Kunkel and 
Chen, 2010]; using radio emissions from the CME core [Tun 
and Vourlidas, 2013]; and using the shock stand-off distance 
from remote observations which has recently shown the pos¬ 
sibility of estimating the field strength upstream of a CME 
[Savani et al, 2011b; Poomvises et ah, 2012; Savani et ah, 
2012]. 

Considering the final focus of estimating the magnetic 
vectors is to predict the terrestrial effects with quantifi¬ 
able uncertainty, the uncertainty in the predicted Kp index 
was estimated by varying the field strength over the range 
I B 1= 18.Ol^^, where cr = 6.9nT [Lepping et al, 2006]. 
The uncertainty in field strength represents the statistical 
average from 82 flux rope fittings estimated between 1995- 
2003. The magnetic vectors were recalculated for each field 
strength and used to drive estimates of Kp which are de¬ 
scribed in more detail in Paper II. 

6. Results 

In order to test the validity of this proof-of-concept archi¬ 
tecture for estimating the magnetic vectors within CMEs, a 
total of eight CME events have been investigated. By using 
the same technique as Eigure 1, the solar disc in 171A AIA 
and projected CME ‘volume of influence’ for these events be¬ 
tween 2010 and 2014 are displayed Eigure 6. Their predicted 
magnetic vectors are displayed in Eigure 7, along with the 
measured LI in situ data. Spherical coordinates are used 
to display the magnetic field rotation over the Cartesian 
system as the orientation components remain independent 
of the magnetic field strength component. The angular ro¬ 
tation and Be) in the predicted magnetic field closely 
follows the broad rotational structure seen within the in situ 
data, with a negative Be indicating a southward magnetic 
field excursion. The deviations in results between the esti¬ 
mated and measured values are discussed below. 

Eor the events investigated, it has been noticed that if 
the overlying magnetic field arcade displayed in solar im¬ 
agery (e.g., within 171A AIA) traverses two active regions 
in close proximity, an adjustment to the standard scheme is 
required. In particular, if the solar arcade is between two 
active regions, the leading polarity is reversed and the ini¬ 
tial magnetic structure is defined by the Bothmer-Schwenn 
scheme from the previous solar cycle. The scenario of this 
more complex behavior is shown in the case study event 
of this article (January 2014, panel G in Eigures 6 and 7), 
as well as in an event on September 2012 (Eigure 6 and 

7, panel E). Therefore, we suggest that the ubiquitous use 
of the Bothmer-Schwenn scheme with a simplistic flux rope 
model is capable of generating a zeroth-order characteriza¬ 
tion of the rotating magnetic field topology with a flux rope 
CME. 

Eigure 7 also illustrates the limitations of a symmetrical 
flux rope model, which is frequently highlighted by a vari¬ 
ety of in situ models, in that the model field strength are by 
definition stronger near the center of the flux rope whereas 
the observed fields occasionally deviate from this pattern. 
As an example, panels (a), (e) and (f) display the strongest 
field near the CME leading edge or sheath, which sometimes 
occurs when a fast CME compresses against the solar wind 
ahead. 

The rotating nature of the magnetic field’s southward ex¬ 
cursion has important consequences for improving start time 


predictions of significant values in Kp index or aid strength 
estimates of the Dst storm onset. Panels (a), (b) of Eigure 7 
show examples of an initial prolonged northward magnetic 
field component and thereby would predict a delayed start 
of large Kp values (see Paper II for more details). 

There are also processes that can influence the accuracy 
of the predicted fields, in particular the interaction of CMEs 
during passage from the Sun to the Earth [Shen et al., 2012]. 
As an example, two CME events launched in quick succes¬ 
sion were detected as a single strong event within in situ 
data and displayed in panel C of Eigure 7. In the interim, 
an experienced observer may be able to manually adjust the 
computational models in response to such unusual situations 
in a heuristic manner used by forecasters. 

The estimation of the impact parameter (i.e., perpendicu¬ 
lar to the flux rope axis) is an important variable in influenc¬ 
ing the predicted magnetic vector. This parameter affects 
the estimated peak magnetic field strength as well the ex¬ 
pected angular change in the field rotation. The influence on 
total field rotation goes from observing a maximum 180° ro¬ 
tation to a minimum of 0° between a trajectory through the 
core and edge, respectively. As an example, the predicted 
vectors for panel C and H in Eigure 7 display a significantly 
larger variation in field rotation than was observed. 

Eor the case of the January 2014 event, draping of the sur¬ 
rounding solar wind magnetic field is likely to account for 
significant portion of the measured terrestrial disturbance 
due to Earth’s trajectory through the outer northern edge. 
As a first principle, the field rotation from a draped mag¬ 
netic field as measured from a 1-D spacecraft trajectory can 
be modeled with the minimal rotations created from a large 
impact parameter modeled flux rope described below. Eor 
cases as extreme as this, a forecast system that generates 
a subtle field rotation may be considered more appropriate 
than generating a ‘missing-Earth’ scenario, but extensive 
statistical analysis will be required to minimize uncertainty 
for such cases. 

Statistically, a spacecraft should have no relationship with 
the CME trajectory, and the frequency distribution of CME 
versus the spacecraft distance is expected to be approxi¬ 
mately uniform. This has not always been observed with in 
situ detectors [Lepping and Wu, 2010], but this is likely due 
to the trajectory being outside their core flux rope behaviour 
[Demoulin et al., 2013]. Therefore a split behaviour of the 
impact parameter, used in this work, between the central 
core and the outer regions is an appropriate choice. A com¬ 
mon uncertainty in impact parameter from various models is 
considered as approximately ±10% [Al-Haddad et al, 2012]. 
Therefore in Paper II, changes to the impact parameter over 
the uncertainty range are used to create an ensemble of pre¬ 
dicted vectors in order to investigate their consequences on 
the Kp index. 

The estimated magnetic vectors from the CME is quasi¬ 
invariant to any trajectory variations that are parallel to 
the flux rope axis. This is because the simplistic model is an 
axis-symmetric cylinder. The small changes to the estimated 
vectors that does occur is a result of small adjustments to 
the radial component of the CME axis direction. The esti¬ 
mated vectors change rapidly once the predicted trajectory 
approaches the legs of the flux rope axis as the influence of 
the radial component is highly non-linear. Under such situa¬ 
tions, detection of a CME with in situ data usually becomes 
inconclusive [Owens et al, 2012] and therefore unlikely to 
have a major impact for the purposes of predicting large Kp 
values at Earth. 

7. Discussion, Conclusions and Future 
Work 

This article displays a reliable mechanism by which mag¬ 
netic vectors can be forecasted. The current process lays 
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the organizational structure that is based on remote sens¬ 
ing and empirical relationships. The example January 2014 
event is severely deflected away from the Sun-Earth line and 
thus highlights the importance of including evolutionary es¬ 
timates of CMEs from remote sensing when attempting to 
provide reliable forecasts (as previously suggested by McAl¬ 
lister et al [2001]). Also, to improve the reliability of the 
magnetic vector forecast, the initial topological structure de¬ 
termined by the Bothmer-Schwenn scheme must be adjusted 
for cases where the overlying field arcade clearly traverses 
two active regions. 

While the current process lays the organizational struc¬ 
ture, in its current format, the concept has not yet been 
statistically proven to be more beneficial at helping esti¬ 
mate the geo-effectiveness of an Earth-directed CME. Eor 
this. Paper II describes a first approach to create an ensem¬ 
ble of magnetic vector predictions that are used to create 
predicted geomagnetic indices, (e.g., Kp). This approach 
leads to a time varying Kp prediction and for those esti¬ 
mates to have quantifiable uncertainties. 

In order to create this proof-of-concept, several assump¬ 
tions and simplifications have been made. This is both a 
strength for the technique being computationally fast, as 
well as a weakness for the simplifications being unable to al¬ 
ways capture the detailed nature of a complicated geomag¬ 
netic storm. A detailed statistical investigation is therefore 
required to further understand the probability distribution 
of accurate forecasts versus false positives. 

The compressed solar wind plasma in between supersonic 
magnetic flux rope obstacles and their driven shock fronts 
has not been addressed in this article, even though they 
have been shown to be significant drivers of magnetospheric 
storms [Huttunen and Koskinen, 2004]. Panel C, D and E 
in Eigure 7 displays the more pronounced examples of high 
amplitude fluctuations in the magnetic field just prior to the 
start of the flux rope CME. Therefore, future work may con¬ 
sider approaches that can better forecast these components 
of geo-effective CMEs. 

In order to determine the usefulness of predicting the 
magnetic vectors for the purposes of estimating geomagnetic 
indices, a standardized procedure that all future techniques 
can be tested against will be beneficial. Such a forecast skill 
score (e.g., Heidke or Brier skill score) will perhaps be more 
useful than a traditional RMSE of individual data points 
between predicted and observed [e.g.. Mays et al.^ 2015], as 
this will potentially prevent uncertainty in arrival time val¬ 
ues skewing the results. 
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Figure 1. Solar source of the 7th January 2014 eruption. 
The Sun is shown as a 171 (Fe IX/X) image from the AIA 
instrument onboard the SDO spacecraft taken at 20.14 
UT before the propagation of the coronal mass ejection 
(CME). The extent of the axis of the CME magnetic 
structure is indicated by the dashed curve, displaying a 
southward deflection from the flare location. The center 
of the axis is shown with a blue square. The volume of 
influence onto the heliosphere from the CME is shaded, 
suggesting that the Earth only grazed the northern edge 
of the CME. The perpendicular distance of the Earth 
from the CME axis is shown in blue. 
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Figure 2. The evolution of the 07 January 2014 CME 
event from three vantage points. The graduated cylin¬ 
drical shell model to estimate the topological structure 
of the event is shown in red. A transition layer ahead 
of the magnetic structure indicates the distance to the 
shock wave driven ahead. 



Figure 3. An empirical model comparing the CME axis 
to Earth distance with a theoretical magnetic structure 
length. The January 7, 2014 event is shown to graze the 
outer edge of the magnetic structure, with a normalised 
impact parameter of Y0=0.91. 
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Figure 4. An empirical model to estimate the radial 
component of the flux rope axis direction (A). The rel¬ 
ative distance of Earth impact and the flux rope nose is 
shown as the red line. The flux rope axis direction is as¬ 
sumed to be perpedicular to the radial at the nose and 
parallel to the radial (90°) at both footpoints. 
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Figure 5. Magnetic vectors from the LI vantage point 
upstream of Earth for the arrival of the CME. The mag¬ 
netic held from the OMNIWeb dataset is shown in GSE 
components (a-c), and spherical coordinates (d-f). The 
red curves overlaid represent the forecasted magnetic vec¬ 
tors at Earth. 
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Figure 6. Solar source of 8 CME eruptions between 2010 
and 2014. The Sun is shown as a 171 (Ee IX/X) image 
from the AIA instrument onboard the SDO spacecraft 
and in the same format as Eigure 1. 
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Figure 7. Predicted (red) and observed (LI; blue) mag¬ 
netic vectors for 8 CME events, where is the angular 
magnetic held direction out of the Sun-Earth plane. 

























































































